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ABSTRACT 

Characterizing the local space density of double degenerate binary systems is 
a complementary approach to broad sky surveys of double degenerates to deter¬ 
mine the expected rates of white dwarf binary mergers, in particular those that 
may evolve into other observable phenomena such as extreme helium stars, Am 
CVn systems, and supernovae la. However, there have been few such systems de¬ 
tected in local space. We report here the discovery that WD 1242—105, a nearby 
bright WD, is a double-line spectroscopic binary consisting of two degenerate 
DA white dwarfs of similar mass and temperature, despite it previously having 
been spectroscopically characterized as a single degenerate. Follow-up photome¬ 
try, spectroscopy, and trigonometric parallax have been obtained in an effort to 
determine the fundamental parameters of each component of this system. The 
binary has a mass ratio of 0.7 and a trigonometric parallax of 25.5 mas, placing 
it at a distance of 39 pc. The system’s total mass is 0.95 M 0 and has an orbital 
period of 2.85 hours, making it the strongest known gravitational wave source 
(log h = —20.78) in the mHz regime. Because of its orbital period and total mass, 
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WD 1242—105 is predicted to merge via gravitational radiation on a timescale 
of 740 Myr, which will most likely not result in a catastrophic explosion. 

Subject headings: Binaries:close-binaries:spectroscopic-stars:individual:WD 1242—105- 
white dwarfs 


Introduction 


As the endpoint of stellar evolution for all stars <8 M 0 , white dwarfs (WDs) provide 
key insight into the late evolution of stellar objects. Double degenerates (DDs), or binary 
white dwarfs, also open a window into the late evolution of binary systems. In particular 
those DDs fated to merge via gravitational wave radiation on timescales shorter than the 
age of the universe represent a potentially significant progenitor population for Type la 
supernovae (SNe la). Closely orbiting DDs additionally r epresen t an impo rtant co mponent 
to the galactic background gravitational wave radiation fjEvans et al.lll987t iHils et al.lll99QT) . 


The first discovery of a close DD system was nearly thirty years ago (iSaffer et ah 


19881 ). with significant efforts to find more via the ESO Super novae Type la Progenitor 


Survey (SPY) and Extremely Low M ass (ELM) WD surveys (jNapiwotzki et al.ll200lL 12003 


Kilic et al.l 120101: Brown et al. 2010). While the SPY survey focused on a large, magnitude 
limited survey for DDs, the ELM survey has focused entirely on apparen tly low mass WDs, 
since these are hard to form without a massive companion (Marsh et ah 1995). Both ap¬ 
proaches have yielded dozens of new DD systems, which allow for statistical investigations 
of such binaries. In this paper, we report the serendipitous discovery that WD 1242—105 is 
a nearby DD system that will merge in a relatively short timescale. 

WD 1242—105 was first reported as a UV-excess sou rc e and misclassifi e d as a subdwarf- 
B star in the Palomar-Green survey flGreen et al.l 1 9861 ). Sa lim &_ Gould (120 02). however, 
identified it as a nearby WD candidate based on its proper motion, which was confirmed 
spectroscopically by IVennes fe Kawkal (120031). with fu r ther observations that did not reveal 
it to be of any particular note (ISalim fe Gouldl l2003t Kawka et al.l 12004]; iKawka &; Vennes 
20061) . Due to its inferred b rightness, gravity, and temperature, it was placed in the loca l 
25 pc sample of white dwarfs (IHolberg et al.ll2008l:lGiammichele et al.ll2012[ISion et al.ll2014l) . 
but lacked any published high resolution optical spectroscopy; presumably it was discovered 


1 This paper includes data gathered with the 6.5 meter Magellan telescopes and the 2.5 meter Dupont 
telescope located at Las Campanas Observatory, Chile. 
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too late for inclusion in the SPY survey. For these reasons, we originally targ eted t his white 
dwarf in a search of photospheric metal line pollution in nearby white dwarfs fjDebes Sz Kilic 


2010 ). The discovery of a clear secondary component to the spectrum around the Ha line 


led us to further investigate the nature of this system. 


In Section [2] we detail the suite of observations we obtained in order to determine the 
nature of this binary. In §3] we analyze our results and place strict constraints on the mass 
and orbital parameters of the binary through simultaneous fitting of synthetic spectroscopy 
and the relative difference in the components’ gravitational redshift. In §5] we place this 
system into the context of other double degenerate systems, and present our conclusions in 

m 


Observations 


2.1. Magellan MIKE spectra 


We o bserve d WD 1242—105 with the blue and red chips of the MIKE spectrograph 
( Bernstein et ah 2003h installed at the 6.5-rn Magellan Clay Telescope at Las Campanas 
Observatory (LCO; Chile) as part of a survey of nearby WDs for the presence of photospheric 
metal lines. All runs used a 0"7x5"0 slit, yielding an average spectral resolution of R ~ 35000 
at 3933 A. The spectra cover wavelengths between 3335 and 9500 A. Each exposure was taken 
with a 600 s integration time to ensure sufficient S/N. Table |T] lists the 27 observations of the 
WD, which were taken over multiple epochs starting in March 2008 with the first discovery 
spectra, and intensive follow-up in April and May 2009. A nearly full period of the orbit 
was obtained in May 2009, allowing us to place precise constraints on the orbital period. 

The data were extracted and flatfielded using the M I KE reduction pipeline written by 
D. Kelson, with methodology described in K elson et ah (120001 ) and K elson (120031 ). Each 
spectrum was corrected for heliocentric motion and each epoch was converted to heliocentric 
Julian date. The continuum around the Ha line was fit with a polynomial and the narrow 
H-a core was used to measure radial velocities for both components of the binary system via 
the simultaneous fitting of two Gaussian curves. In general, this was sufficient to determine 
the radial velocity of the two components at a precision of ~3-5 km/s. When the two 
components of the binary were close to conjunction, the corresponding uncertainty in the 
line centers increased. In those cases uncertainties in fitting each velocity were closer to 
~20 km s —1 . 


The spectrum was also inspected for any evidence of Ca or Mg absorption, indicative of 
accretion due to dust or some external source of metal-rich material. We saw no evidence of 
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this in the raw epoch-to-epoch spectra, but a more detailed analysis is beyond the scope of 
this paper. More stringent upper limits for each component of the binary will be presented 
in a future paper (Todd et al. 2015, in prep.). 


2.2. Time Series Photometry of WD 1242—105 


When the first spectra of WD 1242—105 showed evidence for the presence of a companion 
with a very small orbital separation, we conducted a search for eclipses or any sign of pho¬ 
tometric variability. This included o bservat ions take n on 2009 May 4 using the MagIC-E2V 
instrument on the Baade Telescope fjOsip et al . 2004), a fast readout instrument designed for 
high cadence photometry. In total we took roughly five hours of data, covering nearly 1.75 
orbital periods and using the V filter with exposure times of 60s. The MagIC instrument 
had a field of view of 40"x40" which was large enough to fit both the target and a fainter 
comparison star. A photometric aperture of 25 pixels was used for both WD 1242—105 and 
the comparison. 


2.3. CAPSCAM Astrometry 


We measured the trigonometric parallax using CAPScam OBoss e t al. 2 009! ) with the 
DuPont 2.5m telescope, also at LCO. With a field of view of 6.6', the star and a sufficient 
number of references were observed simultaneously in the standard imaging mode (2048x2048 
pixels, with a pixel scale of 0.194"). Each observing run consisted of taking 15 to 20 exposures 
of 30 to 45 seconds depending on the seeing. Five epochs were obtained between June 
2009 and July 2010 (8th June, 27th January, 10th April, 22nd June, and 31st July). The 
source extraction, source cross matching, g eometric calibration and as t rometric solutio n 
have been obtained using the ATPa software flAnglada-Escude et al.ll2012t iBoss et al.ll2009h . 
The overall precision for the target star and the reference frame stars is 1 mas/epoch. The 
reference stars are selected by the software iteratively based on their epoch to epoch RMS. A 
robust reference frame of 33 stars is used. Parallax and proper motion of the reference stars 
is also obtained as a by-product. Several reference stars have unambiguous USNO B1 and 
2MASS counter parts and the B-K color with the 2MASS magnitudes are used to estimate 
the photometric distance to them. Reference stars with parallaxes <5 mas were used, taken 
either from direct measurements or from photometric estimates. At the end the correction 
from relative to absolute parallax/proper motion is obtained based on 13 reference stars. The 
final parallax and corresponding distance estimation are 25.5±0.9 mas and 39.2 (A 4 pc. The 
statistical uncertainty in the parallax is obtained from by a Monte-Carlo resampling of the 
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astrometry at the same observation dates which properly takes into account all significant 
parameter correlations. The original estimated spectroscopic distance to WD 1242—105 was 
25 pc when it was believed to be a single WD-extrapolating the spectroscopic distance for 
two WDs results in a distance of 35 pc, consistent with our parallax measurement. 


3. Analysis 

In this Section, we describe our analysis of the various observations in order to better 
constrain the properties and fate of the WD 1242—105 system. The combination of all our 
constraints allow for a determination of fundamental parameters for both WDs, which in 
turn allows us to determine whether the system will merge and if so, whether it is massive 
enough to become a progenitor of a Type la supernova. 


3.1. Ephemeris Determination 

Table |T] lists the measured velocities for all epochs of the WD 1242—105 system. The ve¬ 
locities were fit with the non-linear least squares fitting routine curvefit.pro in the Interactive 
Data Language, using sinusoidal curves of the form K cos (pHJD — HJD 0 ) + 7 , with p equal 
to the period, 7 equal to the velocity offset of each component, HJD 0 equal to the reference 
epoch, and K equal to the velocity semi-amplitude. The long baseline of observations as 
well as the dense sampling of the orbit in May 2009 allowed the fitting routines to converge 
on a precise period for the components. These fits give the final ephemeris of the system: 

4 t > = (HJD - HJD 0 )p (1) 

where ifJ£> o =2454970.5901±0.0001, and p=0.118765±0.000002). 

Figure |T| shows the final phased radial velocities for both components, while Table [ 2 ] 
details the velocity semi-amplitudes and 7 for each component of the system as well as 
the inferred mass functions for each component. Because of the two components’ differing 
masses, their velocity offsets represent a combination of the binary’s systemic velocity and 
the gravitational redshift of each component. This can be used in combination with the 
traditional mass-radius relations to solve uniquely for the masses of the two systems, which 
we investigate in §U 
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3.2. Photometric Variability 

Figure [2] shows both the phase folded light-curve of the V band photometry of WD 1242—105, 
with a measured rms of 3 milli-mags and 60-s sampling, along with a Fourier transform (FT) 
of the photometry. No obvious periodicity is seen in the data to a 4-cr level of 1.1 mmag, and 
the overall standard deviation of the photometry matches the estimates of the photometric 
uncertainty. A peak in the FT at 4.28 h is seen, but since this is not coincident with the orbit 
of the system we attribute this to slowly varying atmospheric extinction. The comparison 
star most likely had a different SED than WD 1242—105, and would suffer from differential 
atmospheric extinction which could explain the small variation. Rebinning the data along 
the phase of the orbit does not show any obvious additional structure, implying that the two 
components are well detached, non-eclipsing, and not suffering from any tidal distortion. 


4. Determination of WD properties 

From the orbital radial velocities, spectroscopy, parallax, and photometry we have sev¬ 
eral independent constraints on the mass and radius of each component of the WD 1242—105 
system. We calculate the mass of each component two ways: via the difference in the grav¬ 
itational redshifts of each component, and by simultaneous fitting of the spectra using the 
observed mass ratio and parallax as additional constraints. 


4.1. The mass ratio and gravitational redshift of WD 1242—105 

We can calculate the mass ratio of the two components q = 0.70 ± 0.01 from the ratio 
of the radial velocity semi-amplitudes. The difference in systemic velocity, or 7 provides a 
constraint on the difference of the two objects’ gravitational redshifts: 


_G ( Mi M x 

7 c \Ri qR 2 


( 2 ) 


where Mi and Ri denote the mass and radius of the more massive companion respec¬ 
tively, and i ?2 denotes the radius of the less massive companion. We measure a velocity offset 
difference of 11.6=1=1.3 km s -1 . Under the assumption of theoretical mass-radius relations for 
white dwarfs (which does require an estimate of T e fr), one can un iquel y fit the mas s rat io 
and gravitational red shift to give m a sses o f the two components (jNapiwotzki et al.l l2Q02l) . 
In a recent study by [Holberg et ah (2012), most held DAs showed good agreement with 
theoretical mass-radii relations. Similarly, careful observations of double degenerates show 
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that low mass WDs that have experienced post-common envelo pe ev olution also generally 
behave consistently with expected model mass-radius relations fjBours et al.l 2014). To ob¬ 
tain the mass estimates, we minimized a y 2 metric for the expected Ay for a given primary 
mass and mass ratio, and using the derived T e g- from 1 )4.21 Using this approach, we find 
masses of 0.56^07 and 0.39^0^5 M 0 . We compare this determination of the masses with 
those determined in the next section. 


4.2. Simultaneous Spectroscopic and Photometric spectral energy distribution 

modeling 


Barring a measurement of the inclination of the binary orbit, the unknown masses and 
radii of the binary components need to be disentangled with additional information that can 
be derived from modeling the two components’ optical spectra and their photometric spectral 
energy distributions (SED). This modeling, in concert with the constraints derived from 
the mass ratio and the parallax allows determination of the WD fundamental parameters 
complementary to 1 )4.1) 


In an era of all-sky surveys in the UV through mid-IR, high quality SEDs are now rou¬ 


tine. In particular , WD 1242—105 is within the sky coverage of G ALEX GR7 (IMartin et al. 


2005|), SDSS DR9 ( 

Ahn et al. 

2012 

). 2MASS ( 

Skrutskie et al. 2006 

). and ALLWISE ( 

Wright et al. 


2Qlol ) , resulting in 14 photometric measures of the system’s SED. With an accurate measure 


of the systems parallax, we determine both T e g and log g by fitting synthetic spectra to the 
observed photometry and spectroscopy under the constraints of the observed mass ratio. Fit¬ 
ting the spectroscopic and photometric measurements alone introduces degeneracies where 
multiple similar temperatures and gravities are possible. 


The procedure to fit the spectrum is the same as that used for single WDs where 
the profiles of the hydrogen Balmer lin es are compared to detailed model atmospheres 
(IBergeron et al.lll992l; iLiebert et al.l 120051) . We rely on a combination of spectra taken when 
the components were well separated in velocity space, namely the spectra obtained on HJD 
2454548 (See Table [1]). Each individual (blended) line was normalized to a continuum set 
to unity at a fixed distance from the line center, for both observed and model spectra. 
The atmospheric param eters are t hen found using the nonlinear least-square method of 
Levenberg-Marquardt flPress et ahl 1 9861 ). The uncertainties on fitted values were derived 
from a combination of the covariance matrix of the spectroscopic fitting algorithm, which 
mostly impacts T e ff, and error propagation of the trigonometric parallax and mass ratio un¬ 
certainties, which mostly influence log g values. We determine T e ff for both components from 
the spectroscopic fit but the log g values are fixed from the result of the photometric fit. 

































For our procedure of fitting the photometry, we included the Sloan ugriz photometry 
and converted these measurements into flux densities using the appropriate filters, which are 
then compared with the predic tions from model atmosphere calculations (jBergeron et al. 


1997; Holberg & Bergeron 


20061 ). We apply a correction to the u, i, and z bands of —0.040, 


+0.015, and +0.030, respectively, to acco unt for the offsets between the SDSS filter zeropoints 
and the AB magnitude system (lEisenstein et al. 2006). From the photometry we only fit 
the solid angle with the constraint from the trigonometric parallax. From the observed mass 
ratio and spectroscopic T eff values, we can then determine both gravities assuming a mass- 
radius relation. We iterate on the spectrosopic and photometric fits until we converge to a 
solution on the atmospheric parameters. 


To fit the observations, we re ly directly o n a gr id of mean 3D spectra from pure-hydrogen 
atmosphere 3D simulations ([Tremblay et al. 2013). In this range of T e g, 3D effects on the 
gravities can be quite dramatic. The gravities wer e co nv erted to masse s and radii using 
evolution seque nces with thick hydro gen layers from iFontaine et al.l (120011) for the C/O core 
component and Althau s et al. (120011 ) for the lower mass He core component. The choice of 
composition for the core comes from the implied masses determined in §4.11 He-core WDs 
generally experienced extreme mass loss during the RGB, leaving a core less than 0.5 M 0 . 


Figure [3] shows the resulting model fits to both the optical spectroscopy and the SED 
of the two components under the assumption of the parallax and mass ratio for the system. 
In addition to a good fit to the Sloan photometry, the predicted photometry for GALEX, 
2MASS, and ALLWISE photometry is consistent within the uncertainties. Table [5] lists the 
final T e ff, log g, and masses as derived from our fitting procedure. From this procedure we 
obtain exact agreement in the masses of the two components with that determined in 1 14.11 
This is not completely surprising as the two measurements are linked via the inferred mass 
ratio of the binary system as well as the same theoretical WD mass-radius relationships, but 
does provide confidence in the resulting answer. For our further discussions, we adopt the 
masses and uncertainties derived from the spectroscopic and photometric fitting. 


In summary, the WD 1242—105 double degenerate binary is composed of a C/O core 
WD with a mass of 0.56 M 0 and a less massive 0.39 M 0 Helium core WD. They orbit each 
other in a period of less than 2.85 hours, with an inclination to the line-of-sight of 45.1°. 
The semi-major axis of their orbit is 1 R 0 , which at a distance of 39.2 pc corresponds to a 
maximum angular separation of 120 //as. 
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5. WD 1242—105: The Past, Present, and Future 


The orbital periods of Post-Common-Envelope-Binaries (PCEBs) containing low-mass 
He-core WDs tend to be shorter t han th e orbital periods of PCEBs containing more massive 
C/O core WDs (Zorotovic et ah 2011b The reason for this trend is that shorter period 
systems interact earlier in their evolution and experience enhanced mass loss, ending up 
as lower mass He-core WDs. With an orbital period of 0.118765 d and a 0.39 M 0 binary 
member, WD 1242—105 follows this trend. 

Recently, many close binary white dwarfs with low mass progenitors have been dis- 
covered along with a significant numbe r of longer period and more massive binary systems 
flNelemans et al.ll2005l: iKilic et al.1120121) . with a large number being future merger products. 
WD 1242—105 represents an interesting case as it is close (d<Cl00 pc), with a relatively 
short period and a mass ratio close to 1, but also with a fairly high total mass to the system 
compared to other disc overies. It is a compl ementary detection to both the SPY survey 
( Napiwotzki et al. 2001 1 and the ELM survey ( Brown et al. 2010 1. 


Figure [4] shows the total system mass and merger time for double white dwarfs in 
the ELM Survey and WD 1242—105. The latter has accurate mass measurements for both 
components. We plot the minimum total system mass (and hence the maximum gravitational 
wave merger time) for the ELM white dwarf sample, unless the orbital inclination is known 
from eclipses or ellipsoidal variations. With a merger time of 737 Myr and total mass of 
0.95 M 0 , WD 1242—105 is one of the nearest, most massive, and quickest merger systems 
known. 

There are other double degenerate systems that are likely to be closer than WD 1242—105, 
but lack measured parallaxes or will not merge within a Hubble time. Another merging WD 


syste m with a period roughly twice as long as WD 1242—105, NLTT 53177 (IKarl et al. 


20031 ). may be closer by a few parsecs, given the inferred spectroscopic distance of its two 
components. WD 1242—105 is a near twin of the compact component to the WD 1704+481 
system (IMaxted et al. 2000). which consists of three white dwarfs, two of which are in an 
orbit with a period of 0.145 d. The mass ratio of this pair is also 0.7, with a similar difference 
in their gravitational redshifts. The spectroscopic distance of the distant third component 
is 40 pc (jCfianninas et al.ll201llh which is similar to WD 124 2 —105’s parallax. Fi nally, there 
are other double degenerate systems within 25 pc of the Sun (IHolberg et ah 2008), but those 
have gravitational wave merger times longer than a Hubble time. 

We can also investigate the eventual fate of the system. Fig ure 5 sh ows W D 1242—105 
compared to other DDs and relative to the stability criteria of Mar s h e t al.l 020041 ). which 
dictates whether objects merge violently with the possibility of detonation or stably through 
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Roche lobe overflow mass transfer. All of the massive merger systems in the ELM Survey 
are found due to the « 0.2M 0 ELM white dwarfs. Hence, the se t end to ha ve ext reme mass 
ratios (q ~ 0.2), which should lead to stable mass transfer (Marsh et ah 2004) AM CVn 
objects. On the other hand, WD 1242—105 has q = 0.7, which will lead to unstable mass 
transfer and a merger. 

Simulations of moderately massive C/O core WDs with thick He layers (the case for 
WD 1242—105 when it merges) sh ows that a cataclysmic exp l osion could occur and re¬ 
sult in a SNe la like phenomena flSim et ah 2012; Shen et ah 2010 ) due to the detona¬ 
tion of the helium-shell or double detonation of both the helium layer and the C/O WD. 
Dan et ahl (2 014 ) performed a large parameter space exploration of the merger products for 
CO+CO and CO+He WDs. They find that if the timescale for triple-a reactions is less 
than the dynamical timescale (r nuc < r dyn ), a helium-shell detonation would occur (see also 
Guillochon et ahl 120101 1. Their simulation involving a 0.4 Mq He-core WD with a 0.55 M & 
CO WD (similar to WD 1242—105) reaches a maximum temperature of log T max = 8.4 K and 
density log pT max = 4.65 g cnr 3 . This model has r nuc r dyn , making detonation unlikely. 
The merger will most likely leave behind an extreme helium star (R Cor Bor) with a mass 
close to 0.9 M 0 ( Saio fe Jeffery 2000 b 


WD 1242—105 is also a signicant source of gravitational waves in the mHz frequency 
range. At a distance of 39 pc and i = 45. 1°, we expect th e g ravitational wave strain at 
Earth log/i = —20.78 at log u (Hz) = —3.71 ( Roelofs et al. 2007 b Unfo rtunately, this places 
WD 1242—105 outside of the expected sensitivity of the eLisa mission flAmaro-Seoane et ah 
2012). However, of known ELM systems, it is the strongest source of gravitational wave 
radiation at mHz frequencies, primarily because of its proximity to Earth. 

Finally, we speculate on the discovery of WD 1242—105 and the prospects for finding 
more systems like it within local space. Roughly 5%, and as many as 13% of white dwarfs 
are in close binaries, if one assumes binomial probabilities based on the detection of 2 short 
period DDs within a sample of 44 flMaxted &; Marsh 1999 b The SPY surve y of DA WDs 
found 39 DDs among 679 observed WDs, implying again a 5.7% frequency (Koester et ah 
20091 ). The local sample_of WDs within 20 pc (~126) has 4 reported instances of unresolved 
DDs flHolberg et ahl 120081 ). implying a frequency of ~3% but no more than 7%. These 
numbers are broadly consistent with each other, however there could be at least three more 
local WDs that are actually undiscovered DDs. Given that the 40 pc WD sample should 
include ~70 or so close DD systems, many local DDs are still unaccounted for, but should be 
apparent with the launch of GAIA-these systems will appear over-luminous for their given 
composite gravities, as WD 1242—105 was. Moderate resolution optical echcllc spectroscopy 
of WDs with 10-20 minute cadences, such as what was obtained for WD 1242—105, are 
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sufficient to detect DDs with short periods even with cooler T e g. A volume limited survey of 
double degenerates would provide tight constraints on the degenerate population of binaries 
that may participate in mergers and cataclysmic explosions. 


6. Conclusion 

We have detected a new nearby merging double white dwarf binary system, WD 1242—105, 
previously believed to be a single WD located within 25 pc from Earth. Onr radial velocity 
measurements, photometry, and astrometry show it to be a pair of white dwarfs at a distance 
of 39.2 pc, with a period of 0.1187 d, and possessing a mass ratio of 0.7. The total mass of the 
system is 0.95 M 0 , and since the two components are hydrogen-rich and of similar luminos¬ 
ity we can determine the difference in their gravitational redshifts and thus their individual 
masses. We also simultaneously fit photometry and spectra of the system to calculate the 
individual masses of the binary a second way, which agrees to within the uncertainties. The 
short orbital period of the binary guarantees that it will merge within 1 Gyr, possibly in 
the form of an under-luminous supernova or extreme helium star and makes it one of the 
strongest known gravitational wave sources in the mHz regime. Regardless of its eventual 
fate, WD 1242—105 represents an interesting example of a merging DD system that is bright 
and close to the Earth. 
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Fig. 1.— Phased radial velocity curve of the two components of the WD 1242—105 DD 
system. The red points are derived from Gaussian fits to the core of the Ha feature for the 
less massive component, while the blue points are for the more massive component. The 
solid lines correspond to the best fit orbit for both components, with residuals plotted in the 
lower panel. 
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Fig. 2.— (left) MagIC-E2V high cadence photometry of WD 1242—105 phased to its orbital 
period. Black points show the raw 60s exposures with roughly 3nnnag precision. Overplotted 
in red points are the phase-binned photometry with 9 minute sampling and an rms of less 
than lmmag. (right) Fast fourier transform of the data, showing no significant peaks. A 
long period peak is marginally significant, but does not match the period of the binary. We 
attribute this signal to slow changes in atmospheric conditions. 
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Fig. 3.— (left) Comparison between synthetic spectra of WD 1242—105 and observed spec¬ 
troscopy of the system. Spectral lines from the primary are blue shifted in the spectrum, 
while the secondary’s spectral lines are redshifted. The Balmer series of hydrogen lines are 
compared to simultaneous model fits of each component of the DD system. Black lines are 
the observed MIKE spectra, while red are the model spectra, (right) Comparison of observed 
UV through Mid-IR photometry of the WD 1242—105 system (black error bars) compared 
to the models (black circles). The parallax and mass ratio were used to simultaneously fit 
both components to the spectroscopy and photometry. 
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Fig. 4.— Gravitatio nal wave merger time vs. tota l system mass for double wh ite dwarfs 
in the ELM Survey ( Brown et al. 2013 : Kilic et al. 2014 : Gianninas et al. 2014) and WD 
1242—105 (red filled circle in online version). For the ELM Survey sample, we plot the 
minimum total system mass (assuming i = 90°) when the orbital inclination is unknown, 
and the correct system mass when the inclination is known either from eclipses or ellipsoidal 
variations. Objects to the left of the dotted line will merge within a Hubble time. 
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Fig. 5.— Plot of mass ratio of the sam e syste ms shown in Figure |3j including WD1242, 
compared to the stability criteria of [Marsh et al. ( 20041) . With masses of 0.39 and 0.56 M 0 
and unstable transfer leading to a merger, the resulting merged remnant will likely be an 
extreme helium star, but could result in a sub-Chandrasekhar supernova. 






Table 1. Measured Radial Velocities of WD 1242—105 


UT Date 

HJD-240000 

K, 

km s _1 

a K, 
km s -1 

k 2 

km s -1 

cr K 2 
km s -1 

2008 Mar 23 

54548.72600 

-77 

5 

187 

7 


54548.73416 

-86 

3 

211 

3 


54548.74133 

-75 

3 

179 

6 

2009 Apr 16 

54937.55058 

-8 

8 

87 

6 


54937.64928 

115 

13 

-79 

11 


54937.71134 

-20 

7 

121 

10 

2009 May 10 

54962.47196 

121 

3 

-76 

3 

2009 May 11 

54962.53042 

-41 

4 

143 

4 


54962.53848 

-17 

3 

58 

3 


54962.70545 

135 

3 

-104 

3 


54962.77412 

17 

11 

104 

9 

2009 May 17 

54969.49811 

3 

6 

92 

6 

2009 May 18 

54969.53122 

-66 

5 

188 

4 


54969.59765 

120 

3 

-74 

4 


54969.67401 

0 

3 

74 

3 


54970.57065 

118 

5 

-17 

4 


54970.57816 

-53 

5 

181 

5 


54970.58567 

-76 

3 

208 

4 


54970.59318 

-79 

3 

212 

4 


54970.60070 

-65 

3 

187 

5 


54970.60821 

-26 

3 

121 

3 


54970.61669 

23 

3 

46 

3 


54970.62419 

80 

19 

0 

9 


54970.63170 

113 

3 

-73 

3 


54970.63921 

149 

3 

-123 

3 


54970.64672 

165 

3 

-148 

3 


54970.65423 

163 

3 

-146 

4 
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Table 2. The WD 1242—105 binary 


Parameter 

Primary 

Uncertainty 

Secondary 

Uncertainty 

K (km s -1 ) 

124 

1.2 

178 

1.4 

7 (km s -1 ) 

41.9 

0.8 

30.3 

1.0 

M (M 0 ) 

0.56 

0.03 

0.39 

0.02 

T eS (K) 

7935 

92 

8434 

36 

log g 

7.94 

0.05 

7.54 

0.05 

T C ooi (Gyr) 

1 


0.6 







